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Abstract — The masses of 68 supermassive black holes (SMBHs) in nearby (2 < 0.15) active galactic nuclei 
(AGNs) detected by the INTEGRAL observatory in the hard X-ray energy band (17-60 keV) outside the 
Galactic plane > 5°) have been estimated. Well-known relations between the SMBH mass and (1) the 
infrared luminosity of the stellar bulge (from 2MASS data) and (2) the characteristics of broad emission 
lines (from RTT-150 data) have been used. A comparison with the more accurate SMBH mass estimates 
obtained by the reverberation-mapping technique and from direct dynamical measurements is also made for 
several objects. The SMBH masses derived from the correlation with the bulge luminosity turn out to be 
systematically higher than the estimates made by other methods. The ratio of the bolometric luminosity to 
the critical Eddington luminosity has been found for all AGNs. It ranges from 1 to 100% for the overwhelming 
majority of objects. 
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f INTRODUCTION 

Measuring the masses of supermassive black holes 
(SMBHs) in galactic nuclei is an important part of the 
various studies of both the SMBHs themselves and the 
interrelationship between the evolution of galaxies and 
the growth of black holes in their central regions. 
The most reliable known method of determining the 
SMBH mass involves measuring the kinematics of 
stars or gas motions in the gravitational field of a 
central black hole. Unfortunately, since the angular 
resolution of telescopes is finite, the applicability of 
this method is so far restricted to nearby galaxies: the 
dynamical SMBH mass esti mates have been . obta. ined 
for about 50 galactic nuclei ( Giiltekin et al.|[2009l ). in 
the overwhelming majority normal (inactive) ones. 

Measuring the SMBH masses in active galactic 
nuclei (AGNs) — Seyfert galaxies and quasars — is 
of particular interest, because an active growth of 
black holes occurs in them. The reverberation-mapping 
technique, applicable to type-1 AGNs (with broad 
emission lines in the optical spectrum), is deemed 
to be fairly reliable. It consists in: (1) measuring 
the size of the broad-line region from the delay 
between the line and continuum fiux variations, (2) 
measuring the gas velocity dispersion in this region 
from the line width, and (3) applying the virial 
theorem to determine the black hole mass from these 



two characteristics. The main problem associated with 
reverberation mapping is the necessity of performing 
multiple spectrophotometric observations of AGNs 
for several weeks or months. Therefore, so far the 
SMBHmasses have been measured by this technique 
only in sev eral tens of Seyfert galaxi es and quasars 



(|Kaspi et al.. 2000: Peterson et aljl2004l ). 



Since the methods described above, which may be 
considered direct ones, are difficult to apply, indirect 
methods are also actively used to estimate the SMBH 
masses. For example, in the case of type-1 AGNs, an 
empirical relation between the SMBH mass and broad- 
line characteristics (width an d luminosity) (see, e.g., 
Vestergaard fc Petersonir2006f ) has been widely used in 



recent years. A physical justifica tion for the existence 
of such a relation was offered bv iDibail(ll977f) .In recent 
years, the SMBH masses in AGNs (of both type 1 
and type 2) also came to be often estimated from 
the correlation with the mass or velocity dispersion 
of the host stellar spheroid — an elliptical galaxy or 
the bulge of a spheri c al gal a xv dPerrares e fc Merrittf 



2OOOI: iTremaine et al 
Gultekin et al.ll2009l ). 
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An increasingly active use of indirect methods 
to determine the SMBH masses in AGNs raises 
the question of how reliable such estimates are. In 
particular, is it justifiable to use the correlation between 
the stellar bulge properties and the central SMBH 
mass established for normal galaxies in the case of 
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active galaxies? Previous studies suggested that the 
ratio of the black hole mass to the bulge mass in 
AGNs w as systematical ly lower than that in normal 
galaxies ( Wandel However, it was shown in 

a number of succeeding papers that AGNs obeyed 
the same correlation between the SMBH mass and 
the velocity disp ersion in the bulge (Mr ff-g^) as did 



normal galaxies (Nelson et al."2004 



Onken et al. 2004 



Woo et al. 2010; Graham et al. 20111), but, at the same 
time, significant differences in the correlation between 
the black holemass and th e bulge luminosity were 



revealed ( Nelson et alJ 12004 . iBentz et all I2009D . The 



ity 



interpretation of these results is complicated by the 
fact that the Mbk estimates used for comparison with 
the bulge characteristics were obtained mainly by the 
reverberation-mapping technique and, consequently, 
were determined to within the normalization factor 
/ in the virial formula (see Eq[5]), which can change 
from object to object by several times. In addition, at 
least in the case of one type of AGN, namely narrow- 
line Seyfert 1 galaxies, the fact that the SMBHs have 
masses considerably lower than those expected from 
the correlations with the bulge parameters established 
for normal galaxie s may be con sidered to have been 
firmly established (iMathur et all 120011: I Wandell [2002 : 
Grupe fc Mathuij l2004l) . As an illustration of the 
existing uncertainty, it can be noted that the sample of 
49 galaxies used by ( Giiltekin et al. 2009[ ) to construct 
the correlation between the bulge and SMBH masses 
contains three Seyfert galaxies with highly accurate 
(based on maser emission) SMBH mass estimates, and 
the measured black hole mass in one of them (the 
Circinus galaxy) turns out to be approximately a factor 
of 30 lower than that expected from the correlation with 
the stellar velocity dispersion! 

In this paper, we applied various known methods 
to determine the SMBH mass and the ratio of 
the current accretion rate to the critical one for a 
representative sample of nearby Seyfert galaxies drawn 
from the INTEGRAL hard X-ray all-sky survey. This 
sample has already been used and is being used for 
various systematic studies of the local population of 
AGNs. This work also allows us to reach independent 
conclusions about the applicability of various methods 
for estimating the SMBH masses in AGNs. 

THE SAMPLE 



We used the cata log of AGNs (|Krivonos et al.l 12007 : 
Sazonov et al.ll2007l ). detected in the hard X-ray energy 
band (17-60 k e V) by the IBIS/ISGRI instrument 
(jUbertini et alJ l2003h on board the INTEGRAL 
observatory ( Winkler et al. 2003h during the first three 
and a half years of the mission, from October 2002 
to June 2006. Here, we use a complete sample of 68 
Seyfert 1 (33 objects) and 2 (35 objects) galaxies 



that does not include blazars and objects located near 
the Galactic plane {\b\ < 5°). All AGNs represent 
the local Universe (the redshift of the most distant 
object is 2 = 0.14). The main advantage of the 
sample is that it was drawn on the basis of hard 
X-ray observations and, therefore, essentially does not 
suffer from ordinary selection effects related to the 
absorption of emission in AGNs or the difficulty of its 
detection against the galaxy's bright background. It 
may be considered to be a representative sample of 
Seyfert galaxies in the sense that it properly reflects 
the proportion of type-1 and type-2 AGNs in the 
local Universe and covers almost the entire luminosity 
function of Seyfert galaxies. This sample has been 
used and is being actively used to systematically 
study such properties of the local population of AGNs 
as the luminosity function an d the photoabsor ption 
column density distribution (Sazonov et al.l 120071 ) . the 
total broadband X-r ay spectrum of the population 
( Sazonov et al.l 12008 ). the ratio of the X-ray and 
infrared luminosit ies for AGNs (using data from the 
Spitzer telescope; ISazonov et al.l 120121) , and the ratio 
of the X-ray and narrow-line luminosities (Sazonov et 
al., in preparation). In the subsequent work, it would 
be interesting to use this sample to investigate how 
different AGN emission characteristics depend on the 
SMBH mass and accretion regime. 

Basic information about the objects being studied 
is presented in Table [T] The morphological types and 
redshifts of the galaxies were taken mainly from the 
HyperLeda and NED databases. We took the types 
of AGNs characterizing their optical properties and 



radio loudness from our previous papers (jSazonov et al. 



2007, 2012), where the references to the original sources 
of information are given. The photometric distances 
to the objects were calculated from their redshifts 
for the following cosmological parameters: Hq=72 km 
s^^ Mpc^^, rim=0.3, and r2A=0.7. For several nearby 
objects, we used more accu rate distances from the 
Near by Galaxies Catalogue (|TuI1v et all 120091: ITuIIvI 
Il988l) . In particular, we used the distances to the 
objects to calculate the AGN luminosities in the 17- 
60 keV energy band from the fluxes measured (with 
an accurac y better than 20%) by the IBIS/ISGRI 
instrument ( Sazonov et al.ll2007h . 



ESTIMATION OF THE SMBH MASSES FROM 
EMPIRIGAL RELATIONS 

We applied both well-known indirect methods 
for determining the SMBH mass described in the 
Introduction to our objects. 
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Table 1. Sample of AGNs from the INTEGRAL survey at |6| > 5° 
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Seyfert galaxy, NLSl — narrow-line Seyfert 1 galaxy, BLRG — broad-line radio galaxy, NLRG — narrow-line radio galaxy, XBONG — X-ray 
bright optically normal galaxy; z is the redshift; D is the photometric distance(i/o = 72 km s~^ Mpc~^, n m=0-3, nA=0-7). For several 
nearby objects, we used accurate distances (marked by the asterisk) from the Nearby Galaxies Catalogue JTully et aU[2009l : lTullvlll988D . 
^-17-60 kev is the 17-60 keV (INTEGRAL/IBIS/ISGRI) AGN luminosity. 
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ESTIMATION FROM THE INFRARED BULGE 
LUMINOSITY (BASED ON 2MASS DATA) 

In the last 10-15 years, the existence of a correlation 
between the SMBH mass and properties of the 
spheroidal component of the host galaxy, i.e., the bulge 
of a spiral galaxy or an entire elliptical galaxy, has 
been firmly established. At the same time, the question 
of precisely which characteristic of the spheroid 
(mass, luminosity, velocity dispersion, etc.) correlates 
more closely with the central black holemass remains 
open, because different constructed correlations are 
characterized by approximately the same dispersion — 
the scatter of individual values about the mean (see 
Marconi fc HuntI l2003t iGrahamI 120071 ; iGiiltekin et"ar 



As has already been said, any such correlations 
were firmly established only for normal (inactive) 
galaxies, while their applicability for AGNs is still open 
to question. 

In this paper, we used the correlation between the 
SMBH mass and the near-infrared l uminosity of the 
bulge/elliptical galaxy ( Grahan3l2007[ ): 



IgMsH = -0.37(±0.04)(M/f b + 24) + 8.29(±0.08) , 

(1) 

where A/bh is the black hole mass in solar masses, 
and Mx.h is the extinction-corrected i^-band absolute 
magnitude of the bulge (or elliptical galaxy) . 

The empirical relation ([T]) was constructed by 
Graham from 21 nearby galaxies with available 
dynamical SMBH mass measurements. The scatter of 
IgMsH in this sample about themean (including the 
measurement errors) is 0.33. 

For our purposes, it is convenient to rewrite 
relation ([1]) via the bulge luminosity Lpc.h (using the 
characteristics of the 2MASS Ks filter: the central 
wavelength is 2.16 fim, the bandwidth is 1.686 x lO^'^ 
Hz, and the spectral flux density corresponding to the 
zero magnitude is 666.7 Jy): 



IgAfBH = 0.925 IgL^^b- 31.23, 



(2) 



where Lx^b is measured in erg/s. 

The MsH-iif.b correlation is essentially equivalent 
to the SMBH mass-bulge mass correlation, because the 
infrared luminosity of galaxies is roughly proportional 
to their stellar mass. We preferred to use precisely this 
relation, mainly because almost for all our objects the 
infrared luminosity of a galaxy can be easily estimated 
from the publicly accessible data of the Two-Micron 
All-Sky Survey (2MASS), which covered 99% of the sky 
in three bands: J(A — 1.24 /im ), II{X = 1.66 fim) n 
KsiX = 2.16 fim). 

Most (62) objects of our sample are present either 
in th e Extended Source Catalog (XSC, 
2000 ) or in the Large Galaxy Atlas (LGA, 



Jarrett et al 



Jarrett et al 



{mx) of the entire galaxy is given. This allowed us to 
determine the total infrared luminosities of the galaxies 
(Lk) corrected for interstellar extinction (Table [2]). 

The active nucleus can contribute noticeably to 
the galaxy's infrared luminosity, because part of the 
emission (mainly the ultraviolet one) from the SMBH 
accretion disk is absorbed and reprocessed in the 
surrounding torus of gas and dust. We attempted to 
estimate this contribution. 



^'agn = 



^K,AGN 



(3) 



using the correlation between the infrared and hard 
X-ray AGN luminosities. Such a correlation was 
firmly established, f or example, in our recent paper 
( Sazonov et al. 20121 ) through observations of the same 
AGN sample from the INTEGRAL survey that is used 
here with the Spitzer infrared telescope. It follows 
from the results of this paper and a number of 
published spectral energy distributi ons averaged over 
representati ve samples of quasars ([Elvis et ajj 1994 : 



Sazonov et al. 2004; Ric hards et al.l 12006: Sha ne et al 



201 ih that the K-band luminosity of the dusty torus is 



J/f.AGN 



O.lii 



7-60 kcV, 



(4) 



where Li7_60 kcV is the AGN luminosity measured in 
the INTEGRAL energy band (Table H]). In this case, 
the ratio i/f.AGN /^ly-eo kcV keV can vary from object 
to object by a factor of ~ 2. 

The values of Pagn found from Eqs. ^ and ^ are 
given in Table 2. For the Seyfert galaxy NGC 1068, the 
dusty torus is optically thick for Compton scattering. 
Therefore, the luminosity measured by INTEGRAL 
accounts for only a small fraction (of t he order of 
one percent; see, e.g.. ISazonov et ah 2012 ) of the true 
hard X-ray luminosity of the hot accretion disk corona. 
We estimated this luminosity from Eqs. (15) and (16) 
given below and only then did we calculate Pagn from 
Eqs. dSD and (H. 

The effective angular size of the galaxy in the 2MASS 
image (i?eff)) the radius within which half of the 
galaxy's total K-band luminosity is contained, can serve 
as an additional indicator of the contribution from the 
active nucleus to the galaxy's infrared luminosity. The 
values of i?ofr from the XSC and LGA are listed in 
Table 2. They should be compared w ith the 2MASS 
angu lar resolution, which is ^ 2 ".5 (ISkrutskie et al 



20031 ) of 2MASS, where the apparent K magnitude 



2006h : if Reft is only a few arcsec, then this may indicate 
that the contribution from the active nucleus to Lk is 
dominant. However, in this case, it should be kept in 
mind that the angular size of the galaxy can be small 
per se for the most distant objects of our sample. 

For 17 objects, the contribution from the active 
nucleus to the galaxy's infrared luminosity estimated 
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from the X-ray luminosity turned out to be greater than 
50% (Pagn > 0.5). Therefore, (given some uncertainty 
associated with Eq. @) it turns out to be impossible 
to reliably subtract the AGN contribution from the 
galaxy's total luminosity for such objects, and the 
values of Lx measured for them can be used only 
to obtain upper limits on the SMBH mass from its 
correlation with the bulge luminosity. In addition, in 
all these cases, except for MRK 348 and CYG A, the 
galaxy's size Rcs ^ 5", which also indirectly points to 
a substantial contribution from the central dusty torus 
to Lk- 

The next step in our calculations is to estimate the 
bulge fraction (Pb) m the galaxy's K-band luminosity 
(i.e., basically the bulge mass fraction). Unfortunately, 
an accurate allowance cannot be made for most of the 
sample objects, because many of the galaxies are at 
a distance of the order of or greater than 100 Mpc 
and there are no optical or infrared (2MASS) images 
for them with an angular resolution high enough for 
the nucleus, bulge, and disk emission components to 
be reliably separated. Therefore, we decided to resort 
to a less accurate method, more specifically, to assign 
some mean expected values of Pb to the galaxies 
based on their morphological types listed in Table [T] 
This approach is justified, because the bulge mass 
fraction correlates with the type of galaxy, although 
there exists an appreciable scatter of Pb values among 
galaxies of the same type. There is no uncertainty 
only for elliptical galaxies, for which, obviously, Pb = 
1. For the remaining types, we ado pted the following 
values based on the results from 'Laurikainen et al 



20071 : iGraham fc Wor lev 2008; Laurikainen ct al.l l2010l : 
Pb = 0.25 (SO-Sab), 0.2 (Sb), 0.13 (Sbc), 0.08 (Sc), 
0.06 (Scd and later types). For four galaxies (IGR 
J09446-2636, IGR J10386-4947, IGR J16482-3036, IGR 
J18559+1535) whose morphological types were not 
established but which, nevertheless, are present in the 
2MASS XSC, we adopted Pb = 0.25 as a rough estimate 
of the bulge fraction. 

Thus, we calculated the K-band bulge luminosity for 
most of the objects as 



LK,b — Pb(l - Pagn)Lk, 



(5) 



and obtained the following upper limits for the 
objects in which the AGN contribution is dominant 
(Pagn > 0.5): 

LK,h < PhLii. (6) 

Five objects (IRAS 05078+1626, IRAS 09149-6206, 
IC 4518A, IGR J16558-5203, IGR J17418-1212) enter 
neither into the Extended Source Catalog nor into the 
Large Galaxy Atlas, but they are co ntained only in 
the 2 MASS Point Source Catalog fPSC. ISkrutskie et al. 
20061) . Since the method of calculation described above 



cannot be used for these objects, we made rougher 
estimates of the bulge luminosity for them as follows. 
The PSC provides the apparent K magnitudes of the 
sources measured in an aperture with a radius of 4". 
This angular size corresponds to a linear size from ~ 1.3 
to ~ 5 kpc for the galaxies listed above, i.e., it roughly 
corresponds to the size of the bulges in typical galaxies. 
Therefore, we assumed that the measured magnitude of 
the point infrared source roughly corresponded to the 
luminosity of the galactic bulge, from which, however, 
the expected AGN fraction should be subtracted. Thus, 
we used modified formulas for these five objects: 



PAcr 



L 



_ff,AGN 



LK,h — (1 ^ Pagn)Lk,psc, 



(7) 



(8) 



where Lk,psc is the luminosity of the PSC point source. 

Since the Seyfert galaxy LEDA 138501 is close to the 
Galactic plane in the sky, there are no reliable 2MASS 
photometric data for it. In addition, the morphological 
type of the galaxy is unknown. As a result, this is the 
only object for which even a reliable upper limit on the 
SMBH mass cannot be established from its correlation 
with the stellar bulge luminosity. The SMBH mass 
estimates (or upper limits) for the remaining 67 objects 
are presented in Table [21 

Apart from the A/bh estimates themselves, their 
uncertainties are also important. Since all of the 
objects being studied are bright 2MASS sources, the 
errors related to the photometric measurements of 
their infrared luminosities may be neglected. The 
main uncertainty in our calculations is associated 
with the determination of the bulge fraction in the 
galaxy's infrared luminosity. Based on the results of 
the statist ical studies of galaxies mentioned above 



(see, e.g., Laurikainen et al. 2010l ). we can estimate 



the systematic error in Ig Pb and, consequently, Ig Lx.h 
from Eq. for galaxies with known morphological 
types to be 0.3, i.e., the bulge luminosity can 
be determined to within a factor of 2. For galaxies 
with unknown morphologies, the error in Ig Lx.h is, 
obviously, ~ 0.5— this corresponds to a scatter of Pb 
values from ~ 0.08 to ~ 0.8 about the presumed value 
of 0.25, i.e., almost any type of galaxy is admitted. The 
same large error (~ 0.5) can also be assigned to the 
objects for which the bulge luminosity was estimated 
from the luminosity of the PSC point infrared source. 
Finally, we can assign an uncertainty of ^ 0.1 to the 
values of Ig Lx.h obtained for elliptical galaxies to take 
into account the photometric measurement errors in 
2MASS. Because of the approximate proportionality of 
relation ^ between LK,h and A^bh, the ultimate errors 
in IgAfBH found from the correlation with the bulge 
luminosity are expected to be the same as the errors in 
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Ig Lji^h given above. These uncertainties are taken into 
account in the subsequent statistical analysis. 

Of course, the considerable uncertainty can also 
be associated with the very use of the A/BH^-^K,b 
correlation to determine the SMBH masses in the AGNs 
of our sample, but we can attempt to estimate its value 
only by comparing the mass estimates obtained by this 
method with those obtained by other methods (see 
below). 

ESTIMATION FROM BROAD-LINE 
PARAMETERS (BASED ON RTT-150 
OBSERVATIONS) 

The broad emission lines present in the spectra of 
type-1 AGNs originate in the clouds of photoionized gas 
located inside the sphere of gravitational influence of 
the SMBH. This provides the fundamental possibility of 
using the virial theorem to determine the SMBH mass: 



BH 



/ 



Rv^ 
~G 



(9) 



where R is the characteristic distance from the SMBH 
at which the emission in a particular broad line is 
generated, v is the gas velocity dispersion in this region, 
G is the gravitational constant, and / is a coefficient 
dependent on the structure and orientation of the 
broad-line region. 

The characteristic velocity v is easy to estimate 
from the Doppler broadening of spectral lines. As 
regards the size of the region i?, it can be measured 
by the reverberation-mapping technique. However, this 
suggests a spectroscopic monitoring of AGNs for many 
nights, which is by no means always possible. Therefore, 
in recent years, an indirect method based on the 
estimation of the size of the broadline region from 
the AGN luminosity (L) either in the broadest line 
or in the continuum has been actively used. Dibai 
( I977I ) predicted a power law i? cx L" with a slope 
a = 1/3. It follo ws from reverberation measurements 
that g r-. 0.5-0.7 (IWandel et a l. 1999: Kaspi et aT]l2000l . 
2OO5I : l\^stergaard fc Peterson! i2006, ). The fundamental 
advantage of this method is that the SMBH mass in 
AGNs can be estimated from a single measurement of 
the line width and the line or continuum luminosity. 

We applied the empirical relation between Mbh and 
param eters of the broad Balmer Hg line (A = 4861 A) 
from ( Vestergaard fc Peterson! 2006h to estimate the 
SMBH masses in the type-1 AGNs from our sample: 



Ig Mbh = Ig 



/ FWHM(H^) y 
VlOOO km s'V ^ 



V1042 erg s-i 



+ 6.67, (10) 



where the black hole mass Mbh is in solar masses, 
FWHM is the full width at half maximum of the line, 
and L is the line luminosity corrected for the absorption 
on the line of sight. Relation ([TUl) is characterized by an 
intrinsi c scatter (of Ig Mbh v alues about the mean) of 
- 0.43 (|Vestergaard fc Peter son 2006). 

It is important to note that the virial formula 
([ni is defined to within the scale factor /, because 
there exists a fundamental uncertainty in how the gas 
velocity field is structured in the broad-line regions. 
For exampl e, / = 3 for a spherically symmetric broad- 
line region (IWandel et al. 1999; Kaspi et al. 2000) and 
/ — !■ cx) for a Keplerian disk whose plane is oriented 
toward the observer. Therefore, in several papers, the 
mean (/) for the sample was found in such a way 
that the SMBH masses measured by the reverberation 
technique satisfied the Mbh — cr* correlation for normal 
galaxies. Onken et al. (2004) obtained (/) = 5.5 
precisely in this way. This v alue was then used by 
Vestergaard fc PetersonI I2OO6I to derive the empirical 
formula (fTU)) . Thus, this formula implies that the 
standard Mbh— cr* relation must hold for AGNs. 

In 2008-2010, a series of spectroscopic observations 
was performed with the Russian- Turkish 1.5-m 
telescope (RTT-150) for a sample of 19 Seyfert 1 
galaxies from the INTEGRAL survey with declinations 
5 > —20°. The observations were carried out with the 
TFOSC low- and medium-resolution spectrometer. At 
the first stage of these observations, broadband (3230- 
9120 A) spectra were taken for all objects with a 
resolution of 12 A using grism ^^15 and a 67-/Ltm (1.74- 
arcsec) slit. Subsequently, higher resolution spectra 
were taken for some of the objects in the 3900-6800 A 
band with a resolution of 4 A (grism ^7, a 54 fim 
slit, 1.40 arcsec) and the 5850-8270 A band with a 
resolution of 3 A (grism #8, a 54-/xm slit, 1.40 arcsec). 

We reduced the optical observations and calibrated 
the spectra using the IRAI^ standard software. The 
bias image was subtracted from the two-dimensional 
spectra; subsequently, the images were aligned using 
the flat-field spectra obtained from the spectrum of 
a halogen lamp. The one-dimensional spectra of the 
objects were extracted in a linear aperture that was 
centered on the object's emission maximum and had a 
size large enough for only a negligible fraction of the 
flux in broad Balmer lines to fall outside this aperture. 
The spectral density of the emission was calibrated in a 
standard way using observations of spectrophotometric 
standard stars. 

In most cases, because of the atmospheric jitter, the 
angular resolution during our observations was such 
that the size of a point source exceeded the slit size. 
Therefore, we corrected the fluxes in all of our spectra 
for the slit by assuming the slit to cut off part of the 



^ http://iraf.noao.edu 
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Table 2. Estimates of the SMBH masses from their correlation with the bulge luminosity 



Object 




Rait 






K 


Ak 


IgL 


K 


f AGN 


F 


'b 


Ig 




IgMBH 






(arcscc) 










(erg 


/s) 










(erg/s) 


(Me) 


MRK348 


1 


9 


.0 


10, 


.10 


0. 


.024 


42. 


.70 


0. 


,72 


0. 


,25 


< 


42 


.10 


< 


7.71 


MCG-01-05-047 


2 


23 


.9 


9, 


.39 


0. 


.010 


43 


.10 


0. 


,08 


0. 


,08 




41 


.97 




7.58 


NGC788 


3 


14 


.0 


9, 


.07 


0. 


.010 


43, 


.02 


0. 


,18 


0. 


,25 




42 


.33 




7.92 


LEDA138501 


4 




































MRK1040 


5 


17 


.8 


9, 


.27 


0. 


.035 


43, 


.13 


0. 


,22 


0. 


13 




42 


. 14 




7.75 


IGRJ02343+3229 


6 


24 


. 1 


8, 


.77 


0. 


.036 


43, 


.31 


0. 


,11 


0. 


2 




42 


.56 




8.13 


NGC1068 


7 


15 


.2 


5, 


.79 


0. 


.012 


43, 


.13 


0. 


,30 


0. 


2 




42 


.28 




7.88 


NGC1142 


8 


27 


.0 


9, 


.02 


0. 


.026 


43, 


.71 


0. 


,16 


0. 


2 




42 


.94 




8.48 


1H0323+342 


9 


2 


.4 


11 


.79 


0. 


.078 


43, 


.30 


1. 


,19 


0. 


2 


< 


42 


.60 


< 


8.17 


NGC1365 


10 


59 


.5 


6, 


.37 


0. 


.007 


43, 


.09 


0. 


,01 


0. 


2 




42 


.38 




7.97 


3C111 


11 


3 


. 1 


11 


.38 


0. 


.605 


43, 


.46 


1. 


,44 


1 




< 


43 


.46 


< 


8.97 


ESO033-G002 


12 


6 


.0 


10, 


.03 


0. 


.053 


42, 


.91 


0. 


,17 


0. 


25 




42 


.23 




7.82 


IRAS05078+1626 


13 






*11 


.61 


0. 


.110 


42, 


.29 


2. 


,10 








42 


.29 




7.88 


ESO005-G004 


14 


27 


.5 


8, 


.14 


0. 


.052 


42, 


.59 


0. 


,04 


0. 


2 




41 


.88 




7.50 


MRK3 


15 


8. 


.6 


8, 


.97 


0. 


.069 


43. 


.08 


0. 


,22 


0. 


,25 




42 


.37 




7.96 


MRK6 


16 


2 


.8 


9, 


.56 


0. 


.050 


43, 


.13 


0. 


,21 


0. 


,25 




42 


.42 




8.01 


IGRJ07563-4137 


17 


5 


.2 


10, 


.13 


0. 


.283 


43, 


.09 


0. 


,10 


0. 


25 




42 


.44 




8.02 


ESO209-G012 


18 


4 


.2 


10, 


.09 


0. 


.095 


43 


.62 


0. 


,15 


0. 


,25 




42 


.94 




8.49 


IRAS09149-6206 


19 






*9, 


.43 


0. 


.067 


44. 


.18 


0. 


,10 






— 


44 


.13 


— 


9.59 


MRKllO 


20 


2 


.2 


11 


.80 


0. 


.005 


42. 


.77 


2. 


,74 


0. 


,25 


< 


42 


.17 


< 


7.77 


IGRJ09446-2636 


21 


2 


.5 


12 


.71 


0. 


.031 


43 


.69 


4. 


, 14 


*0. 


,25 




43 


.09 




8.62 


NGC2992 


22 


14 


.2 


8 


.60 


0. 


.022 


42. 


.67 





12 


0. 


,25 




42 


.01 




7.62 


MCG-5-23-16 


23 


5 


.2 


9 


.35 


0. 


.040 


42. 


.51 





47 


0. 


,25 




41 


.63 




7.28 


NGC3081 


24 


21 


.6 


8, 


.91 


0. 


.020 


42 


.60 


0. 


,15 


0. 


,25 




41 


.92 




7.55 


ESO263-G013 


25 


17 


.9 


10, 


.53 


0. 


.062 


43. 


.25 


0. 


,30 


0. 


,2 




42 


.40 




7.98 


NGC3227 


26 


38 


.7 


7, 


.64 


0. 


.008 


42. 


.70 


0. 


,09 


0. 


,25 




42 


.06 




7.67 


NGC3281 


27 


32 


.6 


8, 


.31 


0. 


.035 


43 


.13 





,07 


0. 


,25 




42 


.49 




8.07 


IGRJ10386-4947 


28 


2 


.2 


11 


.82 


0. 


.179 


43. 


.31 





,60 


*0. 


,25 




42 


.71 




8.27 


IGRJ10404-4625 


29 


5 


.4 


10 


.38 


0. 


.058 


43 


.01 





,25 


0. 


,25 




42 


.28 




7.88 


NGC3783 


30 


11 


.9 


8 


.65 


0. 


.044 


42. 


.86 





,30 


0. 


,25 




42 


.10 




7.71 


IGRJ12026-5349 


31 


2 


.6 


9 


.95 


0. 


.075 


43 


.33 





,19 


0. 


,25 




42 


.64 




8.20 


NGC4151 


32 


16 


.6 


7 


.38 


0. 


.010 


42. 


.80 





,38 


0. 


,25 




41 


.99 




7.61 


MRK50 


33 


3 


.5 


12 


.00 


0. 


.006 


42. 


.33 





,73 


0. 


,25 


< 


41 


.73 


< 


7.36 


NGC4388 


34 


36 


.1 


8, 


.00 


0. 


.012 


42. 


.38 





,26 


0. 


,2 




41 


.56 




7.20 


NGC4395 


35 


48 


.5 


9, 


.98 


0. 


.006 


40. 


.47 





, 14 


0. 


,06 




39 


.18 




5.01 


NGC4507 


36 


11 


.9 


8, 


.87 


0. 


.036 


42. 


.99 





,33 


0. 


,25 




42 


.21 




7.81 


NGC4593 


37 


30 


.2 


7 


.99 


0. 


.009 


43 


.13 





,08 


0. 


,2 




42 


.40 




7.98 


NGC4945 


38 


172 


.7 


4 


.48 


0. 


.065 


42. 


.52 





,01 


0. 


,06 




41 


.30 




6.96 


ESO323-G077 


39 


3 


.0 


8 


.80 


0. 


.037 


43 


.23 





,08 


0. 


,25 




42 


.59 




8.16 


IGRJ13091+1137 


40 


9 


.9 


10 


.18 


0. 


.009 


43 


.12 





,36 


0. 


,25 




42 


.32 




7.91 


IGRJ13149+4422 


41 


6 


.0 


10, 


.82 


0. 


.007 


43 


.20 





,41 


0. 


,25 




42 


.36 




7.95 


CENA 


42 


82 


.6 


3, 


.94 


0. 


.042 


42. 


.68 





,02 


1 






42 


.67 




8.24 


MCG-6-30-15 


43 


3 


.2 


9, 


.58 


0. 


.023 


42. 


.33 





,22 


0. 


,25 




41 


.62 




7.26 


MRK268 


44 


4 


.5 


10 


.88 


0. 


.006 


43 


.25 





,35 


0. 


,2 




42 


.37 




7.95 


IC4329A 


45 


2 


.7 


8 


.80 


0. 


.022 


43. 


.28 





,47 


0. 


,25 




42 


.40 




7.98 


NGC5506 


46 


12 


.4 


8 


.19 


0. 


.022 


42. 


.78 





,22 


0. 


,25 




42 


.07 




7.68 


IGRJ14552-5133 


47 


6 


.7 


10 


.36 


0. 


.243 


42. 


.74 





,14 


0. 


,08 




41 


.58 




7.23 


IGRJ14561-3738 


48 


9 


.6 


9, 


.63 


0. 


.031 


43 


.33 





,09 


0. 


,25 




42 


.69 




8.25 


IC4518A 


49 






*11, 


.73 


0. 


.058 


42. 


.10 


1 


,01 








42 


.10 




7.71 


WKK6092 


50 


12 


.1 


10, 


.24 


0. 


.071 


42. 


.70 





,17 


0. 


,2 




41 


.92 




7.54 


IGRJ16185-5928 


51 


4 


.0 


11 


.14 


0. 


.118 


43 


.07 





,40 


0. 


,08 




41 


.76 




7.39 


ES0137-G34 


52 


23 


.7 


8 


.26 


0. 


. 123 


43 


.05 





,03 


0. 


25 




42 


.44 




8.02 


IGRJ16482-3036 


53 


3 


Q 


11 


.34 


0. 


. 124 


42. 


.90 





,71 


*0. 


,25 


< 


42 


30 


< 


7 89 


NGC6221 


54 


42 


.4 


7 


.12 


0, 


.061 


42. 


.88 





,01 


0. 


,08 




41 


.78 




7.41 


IGRJ16558-5203 


55 






*10 


.88 


0, 


.142 


43. 


.58 


0, 


,52 






< 


43 


.58 


< 


9.07 


NGC6300 


56 


54 


.3 


6 


.93 


0, 


.036 


42 


.68 


0, 


,02 


0. 


2 




41 


.97 




7.59 


IGRJ17418-1212 


57 






*11 


.29 


0, 


.210 


43 


.10 


0, 


,62 






< 


43 


.10 


< 


8.64 


3C390.3 


58 


2 


.1 


11 


.54 


0, 


.026 


43 


.30 


2, 


,24 


1 




< 


43 


.30 


< 


8.82 


IGRJ18559+1535 


59 


4 


.3 


12, 


.53 


0, 


.346 


43 


.40 


1 


,53 


*0. 


,25 


< 


42 


.79 


< 


8.35 


1H1934-063 


60 


5 


.3 


9, 


.67 


0, 


.108 


42. 


.60 


0, 


,11 


1 






42 


.56 




8.13 


NGC6814 


61 


31 


.8 


7, 


.66 


0, 


.067 


42 


.81 


0, 


,05 


0. 


,13 




41 


.90 




7.53 


CYGA 


62 


13 


.4 


10 


.28 


0, 


.140 


43 


.85 


0, 


,59 


1 




< 


43 


.85 


< 


9.33 


IGRJ20286+2544 


63 


8 


.1 


9 


.87 


0, 


.164 


42. 


.80 


0, 


,23 


0. 


,06 




41 


.47 




7.12 


MRK509 


64 


2 


.1 


10, 


.01 


0, 


.021 


43 


.47 


0, 


,49 


0. 


,25 




42 


.58 




8.15 


NGC7172 


65 


19 


.0 


8, 


.32 


0, 


.010 


42. 


.87 


0, 


,11 


0. 


,25 




42 


.21 




7.81 


MR2251-178 


66 


1 


.9 


11 


.13 


0, 


.014 


43 


.58 


1 


,18 


1 




< 


43 


.58 


< 


9.08 


NGC7469 


67 


4, 


.7 


8 


.85 


0, 


.025 


43 


.28 


0, 


,14 


0. 


,25 




42 


.61 




8.18 


MRK926 


68 


5, 


.4 


10, 


.64 


0, 


.015 


43 


.49 


0, 


,57 


0. 


,25 


< 


42 


.89 


< 


8.44 



Note. i?cff is the effective angular size of the galaxy (from the XSC and LGA of 2MASS); niK is the apparent magnitude of the entire 
galaxy (from the XSC and LGA), or (m arked by the asterisk) the point source (from the PSC of 2MASS); Ak is the interstellar extinction 
toward the source jSchlegel et al.||1998f) ; Lk is the extinction-corrected luminosity of the entire galaxy or the PSC point source; Pagn is the 
AGN fraction in LK^ Ph is the bulge mass fraction taken to be 0.25 in the absence of morphological information about the galaxy (marked 
by the asterisk); Lk,\> is the bulge luminosity; Mbh is the SMBH mass. 
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Fig. 1. Example of modeling the RTT-150 spectrum. 
The solid curve indicates the spectrum of the Seyfert 
galaxy NGC 7469 near and Hq taken with a 
resolution ^ 4 A. The model (dashed curve) includes 
the broad and narrow Balmer lines as well as the 
narrow forbidden [OIII] and [Nil]. The continuum was 
described by a polynomial and was subtracted. 



flux from a source in the shape of a two-dimensional 
Gaussian. 

The subsequent scientific analysis consisted in 
subtracting the continuum described by a highdegree 
polynomial and modeling the line spectra near h 
H^. The model of the spectrum near H/3 consisted of 
broad and narrow allowed lines (A = 4861.3 A) and 
two forbidden [OIII] lines (Ai = 4958.9, A2 = 5006.8 A). 
The ratio of the [OIII] doublet line fluxes was fixed: 
F2/F1 = 2.91. The profiles of all lines were assumed 
to be Gaussian. The intensity, broadening (measured 
in km/s), and systemic velocity of all the listed lines 
were free model parameters, but the broadenings and 
velocities of all narrow (H^ n [OIII[) lines were assumed 
to be identical (which is reasonable, because these 
lines originate approximately in the same gas) . For the 
spectrum near Hq. we used a similar model that, in 
addition to the broad and narrow Hq lines, included 
the doublet of forbidden [NII[ hues (Ai = 6548.0, 
A2 = 6583.4 A, F2/F1 = 2.96). The broadenings and 
velocities of the narrow and [NII[ lines were assumed 
to be identical. All of the measured line broadenings 
were then corrected for the spectral resolution of the 
measurements. 

The described model allowed the parameters of the 
lines of interest to us to be measured fairly reliably 



for most of the objects (Table |3]). For the narrow-line 
Seyfert 1 galaxy IH 0323+342, we used a simplified 
model in which the Hq and lines had only broad 
components (nevertheless, the complete model was used 
for two other AGNs of the same type, MRK 110 and 
IH 1934-063), because the narrow components of the 
Balmer lines could not be distinguished. An example of 
modeling the spectrum is shown in Fig. [T] The broad 
lines in the spectra of two objects, MRK 6 and 3C 390.3, 
have a complex shape and are poorly described by a 
Gaussian. Therefore, the parameters measured for them 
cannot be considered reliable and were excluded from 
the subsequent analysis. 

The extinction in the path from the broad-line region 
to the observer was estimated from the ratio of the 
Balmer-line fluxes (Balmer decrement): 



Ay.Baim = 7.21 Ig 



f(H„) 
SoF(H/3)' 



(11) 



where i?o = ^o(IIq)/Fo(H^) is the expected Balmer 
decrement for broad lines in AGNs without extinction, 



which was assumed to be 3.06 (|Dong et al.ll2008l) . 

The value of Ay. Balm measured in this way was 
then compared with the interste llar extinction in ou r 
Galaxy toward the object ^y.gai ijSchlegel et al.lll998r ). 
If Av,Ba\in turned out to exceed ^y,gai by more than 
la (the error of our measurement), then we used 
Av = ^v,Baim for the subsequent calculations, i.e., 
we assumed the broad-line emission to be absorbed 
not only in the interstellar medium of our Galaxy 
but also in the object itself (in the galaxy or its 
nucleus). Otherwise, we used Ay — ^y,gai- The final 
correction of the measured line fluxes for extinctio n 
was made by a standard method (jCardelli et al. I ll989l) : 
^(Ha) = 0.8177Ay, A(H^) = 1.1643^y. The derived 
line luminosities are given in Table [3] 

The last column in Table [3] gives the SMBH masses 
calculated from Eq. ([TU)) . The presented errors take 
into account the uncertainties related to both the 
measurement of the H/j parameters and the extinction 
correction. The latter dominates in the objects in 
which an additional (with respect to the extinction 
in the Galaxy) extinction was revealed by the Balmer 
decrement. The uncertainty associated with Eq. (jlOp 
itself was disregarded. 

COMPARISON OF VARIOUS MASS 
DETERMINATION METHODS 

Table [5] gathers together the various estimates 
of the SMBH masses in 68 AGNs from the 
INTEGRAL survey. For 67 objects, we present the 
estimates (Afbuigc) based on the infrared luminosity 
of the bulge/elliptical galaxy obtained from 2MASS 
data here. For 17 Seyfert 1 galaxies, we present 
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Table 3. Estimates of the SMBH masses from Balmer-line parameters 



OGteKT 


Ni 








^V,Balm 




FWHM(H„) 


lgi(Ho 


) 


FWHM(H;3) 


If 


!i(HA 


) 


Ig Mbh 






(10-" 


erg 










(km/s) 




crg/s 




(km/s) 


(erg/s) 


(Mo) 


LEDA138501 


4 


8.22±0 


20 


3 14ib0 


30 







0.53 


5275±144 


42 


82±0 


01 


5981±530 


42 


4g-J-Q 


04 


8.53±0.08 


MRK1040 


5 


6.13±0 


17 


1 07ib0 




1 


96±0.34 


0.32 


3125±80 


42 


21±0 


11 


4349±497 


41 


72ib0 


17 


7.77±0.14 


1H0323+342 


9 


1.18±0 


17 


0.52±0 


05 







0.71 


1879±99 


42 


23±0 


06 


1634±181 


41 


98±0 


04 


7.08±0.10 


3C111 


11 


1.73±0 


03 


0.12±0 


01 


4 


85±0.32 


5.46 


4757±89 


43 


75±0 


01 


6021±608 


43 


34±0 


04 


9.08±0.09 


IRAS05078+1626 


13 


3.33±0 


10 


0.61±0 


08 


1 


81±0.41 


0.99 


4034±132 


41 


95±0 


13 


4549±608 


41 


47±0 


20 


7.65±0.17 


MRKllO 


20 


7.38±0 


59 


2.21±0 


10 





27±0.29 


0.04 


2546±99 


42 


32±0 


10 


2369±141 


41 


80±0 


14 


7.29±0.10 


NGC3227 


26 


24.10±1 


55 


6.52±0 


86 





59±0.46 


0.08 


3504±200 


41 


28±0 


15 


3704±584 


40 


80±0 


22 


7.05±0.20 


NGC4151 


32 


75.80±1 


98 


20.90±1 


50 





53±0.24 


0.09 


5030±160 


41 


75±0 


08 


5430±469 


41 


26±0 


12 


7.67±0.10 


MRK50 


33 


2.47±0 


11 


0.86±0 


09 







0.05 


4064±212 


41 


48±0 


02 


4192±495 


41 


03±0 


05 


7.30±0.11 


NGC4593 


37 


6.74±0 


33 


2.34±0 


26 







0.08 


3645±200 


41 


13±0 


02 


4260±563 


40 


68±0 


05 


7.10±0.12 


IGRJ17418-1212 


57 


1.19±0 


10 


0.20±0 


07 


2 


04±1.07 


1.90 


6137±469 


42 


18±0 


35 


5586±1733 


41 


68±0 


52 


7.96±0.42 


1H1934-063 


60 


3.80±0 


33 


l.OOiO 


14 





68±0.51 


0.97 


1768±64 


41 


27±0 


04 


1418±200 


40 


83±0 


06 


6.24±0.13 


NGC6814 


61 


6.12±0 


31 


2.04±0 


25 







0.79 


2998±146 


40 


84±0 


02 


4015±589 


40 


47±0 


05 


6.91±0.13 


MRK509 


64 


24.00±0 


81 


7.23±0 


47 





26±0.23 


0.19 


4024±127 


42 


86±0 


08 


4580±372 


42 


36±0 


11 


8.22±0.10 


MR2251-178 


66 


10.80±0 


42 


3.75±0 


49 







0.13 


5564±243 


43 


05±0 


02 


6735±791 


42 


61±0 


06 


8.71±0.11 


NGC7469 


67 


22.20±0 


83 


5.89±0 


54 





65±0.31 


0.23 


3205±118 


42 


31±0 


10 


3148±346 


41 


83±0 


15 


7.56±0.13 


MRK926 


68 


8.07±0 


23 


3.00±0 


28 







0.14 


6709±250 


42 


65±0 


01 


9012±1050 


42 


24±0 


04 


8.73±0.10 



Note. j^(Ha ,Hg) are the e xtinction-uncorrectcd line fluxes; Av^Baim is the extinction from the Balmer-line ratio; A^^g^i is the extinction in 
the Galaxy iSchlegel et"im ri998). FWHM (H^, H^) are the full widths at half maximum of the lines; L(Ha, H^) are the extinction-corrected 
line luminosities (we used the boldfaced value of -Av,gal hjih Av.Baim); ^^bh is the SMBH mass. 



Table 4. Cross-correlation of the SMBH masses found by different methods 





Coirelation 


N 


f3±Afi 


a ± Aa 


eo ± Acq 


P 


Pp 


r 


P,. 




-Mbulgc 


8 


-0.12±0.13 


-0.16±0.16 


~ 


0.26 


0.53 


0.06 


0.89 




-Mbulgo 


7 


1.20±0.31 


0.52±0.35 


0.51±0.06 


0.29 


0.54 


0.87 


0.01 




-Mbuige (without NGC 4395) 


6 


0.11±0.30 


-0.05±0.23 


~ 


-0.14 


0.79 


0.18 


0.73 


Mdy„ 




4 


0.46±0.32 


0.01±0.35 


~ 


0.20 


0.80 


0.72 


0.28 




-Mh3 


6 


0.73±0.42 


-0.05±0.30 


0.34±0.05 


0.26 


0.62 


0.60 


0.21 



Note. is tlie number of objects; p is the Spearman rank correlation coefficient; Pp is the probability that the null hypothesis holds; r is the 
Pearson correlation coefficient; Pj. is the probability that the null hypothesis holds. 



the estimates (Mh^) based on the parameters of 
the broad line that we obtained from RTT- 
150 data. In addition, for ten Seyfert 1 galaxies, 
we pro vide the reverberation-based measurements 
(Mr„r) (IPeterson et all 12004 lOnken fc Petersoi] l2002l 



Peterson et alj|2005l iDennev et al.ll2006l) . Fiiiallv. for 

five objects, we give the masses (Mdyn) measured from 
the kinematics of th e stars or gas near the SMBH 
(|Gultekin et al.ll2009l ). 



The various estimates of the SMBH masses in the 
AGNs of our sample are compared in Figs. [5] and [31 To 
get a quantitative idea of how the various methods agree 
between themselves, we calculated the linear regression: 



y = a + /3x , X = \g 



Ml 



y = h 



Mo 







, (12) 



where Mi and M2 are the SMBH masses estimated 
by two different methods. Normalizing the masses to 
IO^Mq (a typical value for our sample) allows the 
correlation between the regression parameters a and 
/3 to be minimized. 

In our calculations, we used the technique described 
in iTremaine et al I I2OO2I and the FITEXY algorithm 
(|Press et all Il992l) . Vlore specifically, the regression 



parameters a and /3 are found by minimizing the 



function x^' 



£2. 

=1 



^0 



2^2 



(13) 



where Xi and yi are the estimates for the ith object, 
exi and eyi are the corresponding errors, and eo is the 
intrinsic scatter of the correlation. The value of eo is 
determined when divided by the number of degrees 
of freedom becomes equal to one. The errors of a, /3, 
and eo are found by varying the parameters within the 
range defined by the condition — Xmin — 1 

Table 2] presents the results of our cross-correlation 
analysis of the SMBH mass estimate from the bulge 
luminosity (Mbuigo), which was used as the variable 
y in Eq. with other types of estimates (Mh/3, 

Mj-cvi Afdyn). The result of a similar comparison of 
Mj-ev with Mh/3 is also presented. In our calculations, 
we disregarded the upper limits on the mass available 
for several objects. 

We found no statistically significant correlation 
between Mrov and Mjj^. This is not surprising, 
because only six objects were compared. However, 
this comparison confirms the validity of the empirical 
formula (fTU|) for Seyfert galaxies, because no systematic 
bias of the M-^p estimates relative to the M^^v estimates 
was revealed, while the scatter of individual values 






108 ^ 



105 ^ ^ 



10^ 10^ 10' 10" 10" 10>» 
MB„(dyn)/M3 

Fig. 2. Comparison of various SMBH mass estimates. 
The regressions in logarithmic space are indicated 
by the sohd hnes. In the middle panel, the dash- 
dotted line indicates the regression withoutNGC 4395 
(object K235). The dashed line indicates the ratio 1:1. 
TheSeyfert 1, 2, and narrow-hne Seyfert 1 (NLSl) 
galaxies are marked by the squares, triangles, and 
circles, respectively. The ordinal numbers of the objects 
are also given in the tables. 



Fig. 3. Same as Fig. [2] for other correlations 



about the mean (eq ~ 0.34 ± 0.05) in our minisample 
agree s well with that found by (jVestergaard fc Peterson 
2006h when calibrating Ea. (|10p based on reverberation 
observations of a larger sample (28 objects) of Seyfert 
galaxies and quasars. 

As regards the determination of the SMBH mass 
from the bulge luminosity, as can be seen from 
Fig. [21 these estimates turn out to be systematically 
higher than those obtained by other methods for 
-^BH ^ IO^Mq (the number of objects for comparison 
is too small for the corresponding conclusion to 
be formulated with regard to higher masses). Our 
regression analysis revealed no statistically significant 
correlation of Afbuigo with Mnp and Mdyn- A weakly 
significant correlation was found between Aij-cv and 
-^buige, with the dependence found being in agreement 
with a linear one (/? = 1.2 ± 0.3). However, only 
one object, NGC 4395, makes a major contribution 
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to this correlatio n. It is often called a "dwarf Seyfert 
galaxy" (see, e.g., Peterson et alj [2005[ ) because of an 
atypically low SMBH mass (Mbh < 10^ Mq) and a 
low luminosity (its luminosity ismuch lower than that 
of the remaining objects in our sample, see Table [1]). 
If we exclude NGC 4395 from consideration, then the 
correlation of M^cv with Mbuigo (based on six objects) 
ceases to be statistically significant. 

The main conclusion that can be drawn from this 
study is that the masses Mbh estimated from the 
infrared bulge luminosity are in poor agreement and, on 
average, are overestimated relative to the more reliable 
estimates of the SMBH masses in Seyfert galaxies. The 
differences of the correlation between the SMBH mass 
and the bulge luminosity for AGNs and normal galaxies 
have al ready been pointed out in several previous 



., pomte 

papers dNelson et al.ll2004HBentz et al.ll2009[ ). 



THE ACCRETION RATE 

Given the SMBH mass estimate, we can calculate the 
Eddington limit. 



iEdd = 1.3 X 10 



46 



IQSMp 



(14) 



and compare it w ith the bolometric A GN luminosity. In 
our recent paper (|Sazonov et al ] |2012|) . based on the set 
of Spitzer infrared observations and INTEGRAL hard 
X-ray observations for the same sample of objects that 
is investigated here, we showed that the hard X-ray 
luminosity is a good indicator of the bolometric AGN 
luminosity, more specifically. 



Li 



7-60 kcV 



(15) 



Table [S] gives the values of LEdd, ^boi and Ly,oi/ L^dd 
obtained for our objects from these formulas using two 
alternative types of SMBH mass estimation: (1) Mbuigc 
(from the correlation with the infrared luminosity of the 
bulge/elliptical galaxy); (2) the best estimate — if there 
is a dynamical mass estimate (Mdyn) for an object, 
then we use it; otherwise (in order of preference), the 
estimate by reverberation mapping (AI^ov), from the 
parameters of the broad H^ line (Mh/3), and from the 
correlation with the bulge luminosity (Mbuigc) is used. 

Formula ([15]) cannot be used for the Compton-thick 
Seyfert galaxy NGC 1068. Therefore, we estimated 
its bolometric luminosity from the infrar ed luminosity 
(Ig 1/ 1 fl = 43.87) of the dusty torus ([Mason et al, 
20061) based on the relation 



Jbol 



Ll5 

by ISazonov et al 



(16) 



20121 from Spitzer 



derived 

observations of Seyfert galaxies from the INTEGRAL 
survey. 



The ratio Lboi/^Edd may be considered as an 
indicator of the SMBH accretion regime. Figure [3| 
presents the values of iboiZ-Z^Edd calculated both from 
the Mbulge estimates (a) and from the best SMBH 
mass estimates (b). For most of the objects, 0.01 < 
iboi/^Edd < 1- This suggests that the gas is most 
likely accreted onto the SMBH in a radiatively efficient 
regime via a ge ometrically thin, optically thick disk 
([Shakura fc Sun vaev 1973) . 

However, the bolometric luminosity for AGN 
IH 0323+342, a narrow fine Seyfert 1 (NLSl) galaxy, 
exceeds the Eddington one. For two more NLSl objects, 
MRK 110 and IH 1934-063, Lboi/^Edd turns out to 
be also high if the best SMBH mass estimate, Mh^, 
is used. In addition, Mbuige for IH 1934-063 exceeds 
Mh/3 by two orders of magnitude. For the remaining 
two NLSl objects, IGR J14552-5133 and IGR J16185- 
5958, -Lboi/^Edd = 0.03 and 0.14, respectively, but these 
estimates cannot be considered reliable, because they 
were obtained from Afbuigc- On the whole, these results 
confirm that the most rapid SMBH growth occurs at 
the present e poch in narrow-line Seyfert 1 galaxies 
( Mathur 2000h and suggest that this growth is far from 
completion. 

For several objects, Lboi/^Edd < 0.01 (Fig. [![). 
This may imply that the gas is accreted onto the 
SMBH at a low rate and/or with a low electromagnetic 
wave emission efficiency. However, in all these cases, 
the values of iboiZ-^^Edd were obtained from the 
Mbuigc estimates and, therefore, may be grossly 
underestimated . 

CONCLUSIONS 

Our study based on a representative sample of nearby 
AGNs has shown that the estimate of the SMBH 
masses in Seyfert galaxies from the infrared luminosity 
(and, consequently, from the mass) of the stellar bulge 
cannot be considered reliable and, on average, leads to 
overestimated values of Mbh (for ^ lO^M©). This may 
indicate that the black hole growth in an appreciable 
fraction of Seyfert galaxies is still far from completion; 
at least we can say about this with confidence in regard 
to NLSl Seyfert galaxies. However, since here we used 
a rather small sample of AGNs and since our estimates 
have considerable uncertainties (in particular, we used 
rather rough estimates of the bulge fraction in the 
galaxy's total luminosity), this conclusion cannot be 
deemed the final one. In addition, it is not yet clear 
how it is reconciled with the conclusion reached by 
several authors (see the Introduction) that AGNs obey 
the standard Mbh-ct* correlation. The ratio of the 
bolometric luminosity to the critical Eddington one 
ranges from 1 to 100% for the overwhelming majority 
of objects. This suggests that the gas is accreted onto 
the SMBHs in Seyfert galaxies at a high rate and 
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Table 5. Summary table of various SMBH mass estimates 



Object 




Is Ma^,^ 

■& dyn 


Ig Mrev 


Ig Mh/3 


J.VJ. bulge 


— f 

.^Edd, bulge 


ibol 


ibol/ 


-7 

.^Edd , best 


ibol/ 






(Mq) 


(Mq) 


(Mo) 


(Mq) 


(erg/s) 


(erg/s) 


-^Edd bulge 


(crg/s) 


-^Edd best 


MRK348 


1 








< 7.71 


< 45.81 


44.51 


> 0.051 


< 45.81 


> 0.051 


MCG-01-05-047 


2 








7.58 


45.68 


43.97 


0.020 


45.68 


0.020 


NGC788 


3 








7.92 


46.02 


44.23 


0.016 


46.02 


0.016 


LEDA138501 


4 












45.29 




46.63 


0.046 


iVl XXXV J. \J'-±\J 








7 77?. 91 
' • ' '7.63 


7 75 


45 85 


44 42 


038 


45 87 


036 










8 13 


46 23 


44 29 


012 


46 23 


012 


IN Vj*^ lUQo 




f\ 046.95 
"■^^6.92 






7 88 


45 98 


44 57 


039 


45 04 


339 




o 
o 






S /I s 
0.40 


4D. 00 


44 88 


U.U.iU 


AP. ^8 
4D.0O 


u.uzu 




9 






T no^ ■ IS 
' -OSggg 


< 8.17 


< 46.27 


45.32 


> 0.114 


45.18 


1.398 


"NT j^i 1 o c r; 


10 






7.97 


46.07 


43.06 


0.001 


46.07 


0.001 


3C111 


11 






no9. 17 

9.088 gg 


< 8.97 


< 47.07 


45.57 


> 0.032 


47.18 


0.025 


ESO033-G002 


12 






7.82 


45.92 


44.09 


0.015 


45.92 


0.015 


IRAS05078-^1626 


13 






7 82 

7.65^ 4g 


< 7.88 


< 45.98 


44.56 


> 0.039 


45.75 


0.065 


ESO005-G004 


14 






7.50 


45.60 


43.13 


0.003 


45.60 


0.003 


MRK3 


15 








7.96 


46.06 


44.38 


0.021 


46.06 


0.021 


MRK6 


16 








8.01 


46. 11 


44.40 


0.020 


46.11 


0.020 


itjKJUfODo-4io/ 


17 








8.02 


46. 12 


44.03 


0.008 


46. 12 


0.008 




18 








8.49 


46.59 


44.74 


0.014 


46.59 


0.014 




19 








^ 9 59 


^ 47 69 


45 14 


^ 003 


^ 47 69 


^ 003 


iVi rtxv i 1 U 


on 






' •^■'7.ig 


^ 7 77 


/t ^ 87 


/I ^ 1 

40. ID 


^ n 1 Q7 


40 .OU 


U.4DO 




21 




< g g2 


< 4g '^2 


46 26 


> Q 35]^ 


< 4g y2 


> 351 


NGn2Qq2 


22 








7.62 


45.72 


43.71 


0.010 


^ 45.72 


0.010 


MCG-5-23-16 


23 








7.28 


45.38 


44.13 


0.057 


45.38 


0.057 


NGC3081 


24 








7.55 


45.65 


43.72 


0.012 


45.65 


0.012 


ESO263-G013 


25 








7.98 


46.08 


44.67 


0.039 


46.08 


0.039 


NGC3227 


26 


7 18l1? 


7.59(2) 


7 OSl oc 


7.67 


45.77 


43.62 


0.007 


45.28 


0.022 


NGC3281 


27 






8.07 


46.17 


43.93 


0.006 


46.17 


0.006 


IGRJ 10386-4947 


28 








< 8.27 


< 46.37 


45.04 


> 0.047 


< 46.37 


> 0.047 


IGRJ 10404-4625 


29 








~ 7.88 


~ 45.98 


44.37 


~ 0.025 


~ 45.98 


~ 0.025 




30 




' 7.38 V-'-.' 




7.71 


45.81 


44.29 


0.031 


45.57 


0.053 


IGRJ 12026-5349 


31 






8.20 


46.30 


44.57 


0.019 


46.30 


0.019 




32 


7 65^-'' 


'■-'-^6. 94 


'■°'7.57 


7 61 


45 71 


44 32 


041 


45 75 


038 


A/ID IT c^n 


33 






7.3O7 19 


< 7.36 


< 45.46 


44. 14 


> 0.048 


45.40 


0.056 


MP Q Q Q 

iNij04o8o 


34 






7.20 


45.30 


43.74 


0.028 


45.30 


0.028 


NGC4395 


35 




^ fir 

5-56g;S(3) 




5.01 


43.11 


41.55 


0.028 


43.66 


0.008 


NGC4507 


36 






7.81 


45.91 


44.46 


0.036 


45.91 


0.036 


NGC4593 


37 




6.99^:«°(4) 




7.98 


46.08 


43.99 


0.008 


45.09 


0.080 


NGC4945 


38 


c 1 r6.32 


6.96 


45.06 


42.49 


0.003 


44.25 


0.018 


ESO323-G077 


39 






8.16 


46.26 


44.08 


0.007 


46.26 


0.007 


IGRJ13091 + 1137 


40 








7.91 


46.01 


44.63 


0.042 


46.01 


0.042 


IGR.I13149+4422 


41 








7.95 


46.05 


44.76 


0.052 


46.05 


0.052 


CENA 


42 








8.24 


46.34 


42.89 


0.0004 


46.34 


0.0004 


MCG-6-30-15 


43 








7.26 


45.36 


43.61 


0.018 


45.36 


0.018 


MRK268 


44 








7.95 


46.05 


44.74 


0.050 


46.05 


0.050 


IC4329A 


45 




< 7.40(1) 




7.98 


46.08 


44.90 


0.067 


< 45.50 


> 0.254 


NGC55UD 


46 








7.68 


45.78 


44.07 


0.020 


45.78 


0.020 


H-jiXj i4ooz-o ioo 


A T 
4 / 








7 OQ 


45.33 


43.84 


0.033 


45.33 


0.033 




48 








8 25 


46 35 


44 22 


008 


46 35 


008 


IC4518A 


49 








< 7.71 


< 45,81 


44.06 


> 0.018 


< 45.81 


> 0.018 


WKK6092 


50 








~ 7.54 


45.64 


43.89 


0.018 


45.64 


0.018 


IGRJ16185-5928 


51 








7.39 


45.49 


44.62 


0.137 


45.49 


0.137 


ES0137-G34 


52 








8.02 


46.12 


43.43 


0.002 


46.12 


0.002 


IGRJ 16482-3036 


53 








< 7.89 


< 45.99 


44.70 


> 0.052 


< 45.99 


> 0.052 


NGC6221 


54 








7.41 


45.51 


42.89 


0.002 


45.51 


0.002 


IGRJ16558-5203 


55 








< 9.07 


< 47.17 


45.24 


> 0.012 


< 47.17 


> 0.012 


NGC6300 


56 








7.59 


45.69 


43.02 


0.002 


" 45.69 


" 0.002 


IGRJ17418-1212 


57 






7.96?;^^ 


< 8.64 


< 46.74 


44.85 


> 0.013 


46.06 


0.062 


3C390.3 


58 




8.46^:11(1) 


9.61«;3 


< 8.82 


< 46.92 


45.60 


> 0.048 


46.56 


0.111 


IGRJ18559+1535 


59 




< 8.35 


< 46.45 


45.53 


> 0.122 


< 46.45 


> 0.122 


1H1934-063 


60 






6.24|;?I 


8.13 


~ 46.23 


43.58 


^ 0.002 


" 44.34 


^ 0.176 


NGC6814 


61 






6 gr"* 


7.53 


45.63 


43.42 


0.006 


45.01 


0.026 


CYGA 


62 


q ^o9.53 
^■^^9.28 




< 9.33 


< 47.43 


45.57 


> 0.014 


47.53 


0.011 


IGRJ20286+2544 


63 






7.12 


45.22 


44.11 


0.079 


45.22 


0.079 


MRK509 


64 




8.16«;?«(1) 


Q qq8 ■ 32 
"■^^8.12 


8.15 


46.25 


45.11 


0.073 


46.26 


0.072 


NGC7172 


65 




7.81 


45.91 


43.87 


0.009 


45.91 


0.009 


MR2251-178 


66 






8-7i|:f 


< 9.08 


< 47.18 


45.60 


> 0.027 


46.81 


0.062 


NGC7469 


67 




7.09^;J^(1) 


7.56^;6| 


8.18 


46.28 


44.38 


0.013 


45.19 


0.157 


MRK926 


68 




Q "708 . 83 

0. (Og^ea 


< 8.44 


< 46.54 


45.20 


> 0.046 


46.83 


0.024 



Note. Mdy n is t he dynamical mass estimate (Giiltckin ct al."20Q9); A'/^-cv is the rcvcrberation-based mass estimate from: (l)lP eterson et all 
I I2004I') . (2)|Onken & Peterson (2002), (3) Peterso n et al. ( 2005 ), (4) Denney et al. (2003); Mh^ is the mass estimate from the parameters of 
the broad Une (RTT-150); A/buigc is themass estimate from the bulge luminosity(2MASS); Lboi is the bolometric luminosity; -£/Edd, bulge 
is the Eddington luminosity calculated from A/tuigej ^Edd,best is the Eddington luminosity calculated from the best mass estimate. 
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in a radiatively efHcient regime. This conclusion is 
consistent with the resu l ts of previous s t udies (see, e.g., 
Middleton erall (|2008[ ). lde Rosa et all (|2012l )). 



ACKNOWLEDGMENTS 

We wish to thank the referee for a number of 
useful remarks. We used data from the Russian- 
Turkish 1.5-m telescope, the 2MASS survey (University 
of Massachusetts, IPAC), and the NASA/IPAC/NED 
(Jet Propulsion Laboratory, California Institute of 
Technology), VizieR (CDS, Strasbourg), and Hyper- 
Leda (Observatoire de Lyon) databases. This study 
was supported by the Russian Foundation for Basic 
Research (project nos. 09-02-00867-a, 10-02- 01442-a, 
ll-02-12285-ofi-m-2011, 11-02-12271- ofi-m), Programs 
P-21 and OFN-16 of the Russian Academy of Sciences, 
and the Program for Support of Leading Scientific 
Schools of the Russian Federation (NSh-5069.2010.2). 



1 r- 



0.1 =- 



10- 



10- 



10-1 



t.., = 0.01 



bol ^■^^'-'Edd 



10^ 



10' IC 



10« 



10'" 



Fig. 4. Top: Bolometric-to-Eddington AGN luminosity 
ratio versus SMBH mass. AIbh was estimated from 
the infrared bulge luminosity, except for LEDA 138501 
(marked by the four-point star) for which the estimate 
from the parameters of the broad line is presented. 
The Seyfert 1, 2, and narrow-line Seyfert 1 (NLSl) 
galaxies are marked by the squares, triangles, and 
circles, respectively. The measurement uncertainties are 
not shown. Bottom: The same, only based on the 
best mass estimates: Mdyn (six-point stars), M^cv (five- 
point stars), Mjj^ (four-point stars), Mbuigc (remaining 
points). 
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